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Abstract 

Beginning with a set of simplified models for spin-0, spin-^, and spin-1 dark matter candidates 
using completely general Lorentz invariant and renormalizable Lagrangians, we derive the full set of 
non-relativistic operators and nuclear matrix elements relevant for direct detection of dark matter, 
and use these to calculate rates and recoil spectra for scattering on various target nuclei. This 
allows us to explore what high energy physics constraints might be obtainable from direct detection 
experiments, what degeneracies exist, which operators are ubiquitous and which are unlikely or sub¬ 
dominant. We find that there are operators which are common to all spins as well operators which 
are unique to spin-| and spin-1 and elucidate two new operators which have not been previously 
considered. In addition we demonstrate how recoil energy spectra can distinguish fundamental 
microphysics if multiple target nuclei are used. Our work provides a complete roadmap for taking 
generic fundamental dark matter theories and calculating rates in direct detection experiments. 
This provides a useful guide for experimentalists designing experiments and theorists developing 
new dark matter models. 
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I. INTRODUCTION 


The existence of non-baryonic dark matter has been inferred from measurements includ¬ 
ing galactic rotation curves [1], large scale structure surveys [2-4], X-ray observations [5], 
gravitational lensing [6, 7], and cosmic microwave background anisotropy measurements [8], 
spanning cosmological eras from the present day to the remote past. This widespread and 
robust data has led to cold dark matter models with a cosmological constant, labeled ACDM 
becoming entrenched as the standard cosmological model. 

Nevertheless, this impressive array of observations has only been sensitive to the grav¬ 
itational influence of dark matter and constrained its relic abundance, leaving its particle 
nature as one of the most important open questions in physics. The search for dark matter 
includes indirect astrophysical searches ([9-13]), collider production efforts (for some exam¬ 
ples of dark matter searches at the LHC, see [14-18]) which will examine new territory soon 
with LHC run 2 which will commence this year, and attempts to observe dark matter inter¬ 
actions with Standard Model (SM) particles via dark matter-nucleus scattering processes in 
direct detection experiments, to which we now turn. 

The search for dark matter via direct detection goes back at least three decades [19, 20] 
and has been particularly vigorous over the last decade or so with experiments such as LUX 
[21], XenonlOO [22], CDMS H (Ge) [23], CDMS I (Si) [24], DAMA/LIBRA [25], COGENT 
[26], and GRESST [27] pushing ever deeper into weakly interacting dark matter mass and 
scattering cross-section parameter space, but has thus far failed to yield a convincing signal. 
In the near future detectors such as Super GDMS [28] (which has recently released its hrst 
results on low mass dark matter searches [29, 30]), XENONIT [31], and DARWIN [32] are 
expected to push the limits of direct detection orders of magnitude below the current levels. 

In order to connect observations to microphysical models one needs a general framework 
within which to interpret the observations of direct detection experiments. For quite some 
time the prevailing method of analyzing dark matter-nucleus interactions has been to assume 
that dark matter is a weakly interacting massive particle (WIMP), and then to categorize the 
interactions as elastic and isospin conserving and either spin-independent or spin-dependent 
[33, 34]. For some well studied models of dark matter, such as the weakly interacting 
Majorana neutralino found in supersymmetry models, this assumption is reasonable. 

With an absence of observed dark matter signals, there has of late been a surge in 
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interest in exploring more general types of interactions between dark matter and nuclei. 
Generalizations include inelastic and momentum dependent interactions, which may arise 
due to additional structure in the dark sector including excited dark matter states, or dark 
gauge bosons giving rise to electric and magnetic form factors [35-42], 

The formalism of choice for many of these investigations is relativistic effective held 
theory, which provides a model independent framework to analyse dark matter-SM inter¬ 
actions [43-45]. It has been shown that these effective theories break down when applied 
to high-momentum transfer experiments, such as the LHC [46]. Therefore analyses moved 
beyond this framework and have moved to what are labeled as ‘simplihed models’ in¬ 
stead [47-49]. Simplihed models are held theories which extend the SM by a single dark 
matter particle and a single mediator particle which allows the WIMP to communicate with 
quarks and/or leptons. The newly added dark matter content is assumed to be a singlet 
under the SM gauge groups (we will consider some cases where the particles mediating 
the interaction have SM charge). In this context it is then possible to calculate collider 
amplitudes valid at the high energies of interest in such experiments. Given this simple dark 
sector, one can write down an exhaustive list of every combination of WIMP and mediator 
spins, and all possible tree level interactions. These simplihed models have now gained 
popularity for analyzing indirect detection signals [50, 51], allowing connections to be made 
with the growing body of literature which make use of them. 

Another step towards placing dark matter-nucleus interactions on a general footing has 
been accomplished recently by utilizing a non-relativistic ehective held theory (EFT) ap¬ 
proach [52-55]. Since the interactions in direct detection scenarios are assumed to take 
place due to an incoming dark matter particle with a typical velocity (P(100km/s), the 
recoil momenta in such an interaction will be 0{'^ lOOkeV). The particle masses involved, 
including the nucleons of roughly GeV scale, the dark matter particles, which typically range 
from the GeV region to several orders of magnitude above, and mediators that can also be 
quite heavy compared to the typical interaction momenta, produce a situation where an 
EFT treatment is quite natural. 

In order to circumvent as much model dependence as possible, one can construct general 
interactions which obey Galilean invariance, T-symmetry, and Hermiticity. These operators 
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will take the standard effective four-particle interaction form, reminiscent of Fermi’s original 
model of weak interactions. The non-relativistic interactions can be shown to be functions of 
only four parameters including the nucleon spin S^, the dark matter spin S^, the momentum 
transfer, q, and a kinematic variable v~^ which is a function of the relative incoming (v^^in — 
VN,in) and outgoing velocities v^^out - VN,out 


=2 


VN,in T '^x^out '^N,out) Vjx^in T 


q 


2/i7v 


( 1 ) 


which obeys ■ q = 0. It was demonstrated in [53] that there exist fifteen such non- 
relativistic interactions which arise from twenty possible bi-linear combinations of dark 
matter and nucleons. 

The formalism developed in [53] is unique in being the only analysis to comprehensively 
develop the nuclear physics of direct detection experiments. From this general framework it 
is now apparent that there are interactions beyond the standard spin independent/dependent 
type. The origins of these ‘new’ interactions are not necessarily exotic and it has been shown, 
in the context of relativistic EFT, how many of them can be generated [56]. 

What has been lacking to date however, is a completely general and comprehensive treat¬ 
ment that connects high energy microphysics with low-energy effective nuclear matrix ele¬ 
ments in a model independent way. It is possible, for example, that the various interactions 
listed in [53] can give rise to degeneracies where different fundamental dark matter La- 
grangians, describing dark matter and interaction mediators of various spins, can produce 
the same interaction types. This will obviously pose problems for attempts to discern the 
properties of dark matter when interpreting the results of experimental data. Furthermore, 
dark matter may not be spin-^, which creates a need for extending the parametric frame¬ 
work from the four descriptors listed above. In particular, as we shall show, this allows the 
existence of new non-relativistic operators to appear in the low energy effective theory. 

Motivated by the above we present here a general analysis covering a broad spectrum 
of particle and interaction types, starting from the microphysics, which will enable one to 
link experiment with fundamental theory while incorporating the new nuclear responses 
described in [53]. 

In this work we build upon the NR-EFT description by examining simplified models which 
incorporate the most general renormalizable Lagrangians for scalar, spinor, and vector dark 
matter interacting with nucleons via scalar, spinor, and vector mediators, consistent with 
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Lorentz invariance and hermiticity while imposing stability of the dark matter candidates. 
We integrate out the heavy mediator and obtain effective relativistic interaction Lagrangians. 
Next, we take the non-relativistic limit of these Lagrangians, and identify them with the NR 
operators from [53], which are reproduced below, in Table 1. Using these, we identify which 
electroweak nuclear responses are excited by a given fundamental interaction model and 
determine the relative importance of various models within the context of direct detection 
experiments consisting of xenon and germanium targets by exploring the relative magnitude 
of coefficients of these operators, and also their energy dependence. 

The paper is organized as follows; in section II the EFT formalism of [53] is summarized, 
in section III we build the generalized relativistic Lagrangians and in section IV we out¬ 
line the signatures and distinguishability of these models in the context of direct detection 
experiments, providing a framework for both experimentalists and theorists to base their 
future analyses. 


II. EFFECTIVE FIELD THEORY OF DIRECT DETECTION 


Conventionally, coherent WIMP-nucleus scattering has been considered to come from two 
types of interactions; spin-independent (SI) and spin-dependent (SD). SI interactions couple 
to the charge/mass of the nucleus while SD couples to the spin. The nuclear cross section 
is generally written in terms of the nucleon cross section at zero momentum transfer, Uo, 
and a form factor, F{q), to take into account the loss of coherence over the finite size of the 
nucleus. 


da 

dEr 


M 



(g) + . 


( 2 ) 


where M is the mass of the target nucleus and is the WIMP-nucleus reduced mass. 
This picture has recently been shown to be incomplete, as it is also possible for the WIMP 
to couple to the nucleus through additional nuclear responses [53] . Working in the language 
of a non-relativistic (NR) effective field theory Fitzpatrick et ah identified 15 operators 
to characterize the ways in which a WIMP can couple to the various nuclear responses. 
These operators are constructed from combinations of non-relativistic vectors which respect 
Galilean invariance, T symmetry and which are Hermitian. We list them in table 1. The 
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Hermitian vectors are: 


i —, v^ = v + ^, Sn, (3) 

where q = p' — p = k — k' is the momentum transfer, v is the velocity of WIMP with respect 
to the nucleus of the detector, Pat is the reduced mass of the system and and are the 
WIMP and nuclear spins respectively. Throughout the paper, we denote by p and p' the 
incoming and outgoing WIMP momenta and by k and k' the incoming and outgoing nuclear 
momenta respectively. Energy-momentum conservation implies the orthogonality condition 
q-v^ =0. Here we will briefly outline the procedure employed in [53] in going from the NR 
operators to the hnal differential WIMP-nucleus cross section. 


TABLE I. List of NR effective operators described in [53] 


Oi 


02 


03 

X 

04 

S^-Sn 

05 


Oe 


O7 

Sn ■ v-^ 

Os 

Sx-v^ 

O9 

*4 • (Sn X 
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Sx ■ {Sn X v-^) 
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■ Sn) 

Ol4 

^(Sn-v^){^-Sx) 

Ol5 

-(4 ■ A) ((®«X ■ 


g 

mjv , 


In general one can write down the non-relativistic interaction Lagrangian as 

( 4 ) 

a=n,p i=l 
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where the coefficients cf are given by the microphysics of the interaction and in general one 
could allow for isospin violation by having different couplings to neutron and proton inside 
the nucleus. This can be rewritten in 2-coniponent isospin space as 

(5) 

T=0,1 i=l 

where and are the identity matrix and the Pauli matrix respectively. The nucleus is 
composed of nucleons, and these can individually interact with the WIMP. This is incorpo¬ 
rated by considering the operator 0{j) as an interaction between a single nucleon, j, and 
the WIMP, and then summing over the nucleons. 

15 15 A 

z '£c'o.t-^ Y. (e) 

r=0,l i=l T=0,1 i=l j=l 

where A is the atomic mass number given by the total number of neutrons and protons. 
One can do the same reduction with A, 

= - {p7v(i),i = 1,..., A - 1} (7) 


where Vy. and VNi'i) are the symmetrized combination of incoming and outgoing velocities 
for the WIMP and nucleons respectively, (here T stands for target, i.e., the nuclear 
center-of-mass) is dehned as 



1 


A 

'y ( \yA r^7V,ont(^)] 

i=l 


( 8 ) 


This allows for a decomposition of the nucleon velocities into internal velocities hjv(i) that 
act only on intrinsic nuclear coordinates and ‘in’ and ‘out’ velocities that evolve as a WIMP 
scatters off the detector. As an example, the dot product between and Sn can be 
rewritten as 


^ ^ 1 ^ 

^ y y ( — ' Sixil) (9) 

i=l ^ 

^ - r ^ 1 

= py ■ ^ SN{i) “1^2 + VN,out{i)] ■ SN{i) \ (10) 

i=l \ i=\ J 2nt 


The second term in the curly brackets is internal to the nucleus and acts as an operator on 
the ‘in’ and ‘out’ nucleon states. V]v,in can be replaced by pN,in/M acting on the incoming 
state, which can in turn be replaced by /M, and similarly PN,out/M by /M on the 
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outgoing nuclear state. Finally, since the nucleus is non-zero in size and individual nucleons 
locally interact with the WIMP, nuclear operators built from Oi are accompanied by an addi¬ 
tional spatial operator where x{i) is the location of the nucleon inside the nucleus. 


Starting from Eqn. 6 and using the substitution rules for n-*“ and including a factor of 
^-iq-xi^ the interaction Lagrangian can be written as a sum of hve distinct terms (nuclear 
electroweak operators) that only act on internal nucleon states. Their coefficients, on the 
other hand, act on WIMP ‘in’ and ‘out’ states. The WIMP-nucIeus interaction can then be 
written as 




where 


r=0,l 
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The WIMP-nucleus amplitude, M, can then be succinctly written as 


M= (Jx, ^x;jN, Mjvl {/o"^ + + 5 . P + 4 . g + 4 . M^). 

r=0,l 

(14) 

By using spherical decomposition, the internal nuclear operators S, T, P, Q and R can be 


further rewritten in terms of standard nuclear electroweak responses as follows: 



Where there is an implicit sum over the nucleons. 


A 

Ojm-Aq) = (16) 

i=l 

and the various electroweak responses are dehned as 


MjuiqA = jj^qAAni^x) 

^JL = jj{.qx)YjLMiAx) 

Ajm = Mfj{qxi) ■ 

^'jM = |^ViMjA^(gfi)| ■ a{i) 

^'jM = ^Vi X Mfjiqxi) ■ <?(*) X ^ Vi + }^MAj(qXi) ■ a{i) 

_q j [ q j z 

^'jM = i ^ViMjM(gTi) ■ (T(i) X ^Vi 

[q J L ? J 

Eja^ = Mfj{qXi) ■ a{i) 

Qjm = VtjM{qXi) + -T!'j^{qxi) 

^JM = ^jM{qXi) - ^SjA^(ga;i) 

Km = KMi.qXi) - ^MjMiqXi) 


(17) 
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where Yjm and Yjlm are spherical harmonics and vector spherical harmonics respectively. 
We are only considering elastic transitions, and assnming parity and CP as symmetries of the 
nuclear ground state. This eliminates some of the responses, and only M, S', A, S”, 
survive. To calculate cross-sections, one needs to square the amplitude, average over initial 
spins and sum over hnal spins. The matrix element squared for the nuclear portion of the 
amplitude has been made available by Fitzpatrick et al. [53], and codes have been supplied 
to calculate the full amplitude and rate [54]. 

As we shall describe, in the following analysis we discovered that two additional NR 
operators are required to fully describe the scattering of spin-1 WIMPs off nuclei, 

rriN 

0,s = t—-S-SN, (18) 

niN 


where S is the symmetric combination of polarization vectors. Appendix A contains the 
details required to include these new operators in the above formalism. 


III. SIMPLIFIED MODELS FOR DIRECT DETECTION 

From a model building perspective, one would like to know how relevant the novel nuclear 
responses are in interpreting direct detection data. Previous work [56] demonstrated that 
using only the SI/SD form factors (even with additional momentum dependence taken into 
account) can lead one to infer wildly incorrect values of the WIMP mass and cross sections. 

Here we go further by starting with simplified models at the Lagrangian level, where 
‘simplihed model’ means a single WIMP with a single mediator coupling it to the quark 
sector. This is useful for two reasons; it allows us to better explore which NR operators 
arise from a broad set of UV complete theories, and also make connection with the growing 
body of literature which use simplihed models for indirect detection and collider searches. 

When it comes to interpreting signals, knowing comprehensively how different interac¬ 
tions with different nuclei arise from different UV complete models will allow us to identify 
degeneracies between competing models. Further, it can also help optimize target selection 
for maximum discrimination of the UV model parameter space. 
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In building these simplified models we remain agnostic about the WIMP’s spin, and 
consider dark matter spins of 0, | and 1. We do however only consider renormalizable inter¬ 
actions between quarks and WIMPs. To ensure a stable WIMP, we assume that the WIMP 
is either charged under some internal gauge group or a discrete symmetry group (for example 
Z 2 ). However, we assume that this gauge charge is not shared by quarks. We will couple 
the WIMP to the quarks via a heavy mediator in two distinct ways: charged and uncharged 
mediators, each with all possible spins consistent with angular momentum conservation. 
The mediator mass is chosen to be the heaviest scale in the problem (and certainly much 
greater than the momentum exchange which characterizes the scattering process) so that we 
can integrate it out (see appendix B for details). This leads to relativistic effective WIMP- 
nucleon interactions, whose NR limit can then be examined. In the uncharged mediator case 
we will consider mediators that are neutral under all SM and WIMP gauge charges, while 
in the charged case, the mediator must have both WIMP and SM gauge charges. Given the 
above as a guide, our Lagrangian construction is then constrained only by gauge invariance, 
Lorentz invariance, renormalizability and hermiticity. In certain cases which follow, the re¬ 
quirement of hermiticity demands coupling constants be complex. Unless explicitly noted, 
the coupling constants are dimensionless and can be assumed to be real. 


A. Uncharged-mediator Lagrangians 


1. Scalar Dark Matter 


We begin with a spin-0 scalar WIMP, S, which has some internal charge to ensure sta¬ 
bility, and S'! is its Hermitian conjugate. To have renormalizable interactions, the neutral 
mediator can only be a scalar or a vector. We denote the scalar mediator by (j) and the 
vector mediator by with field strength tensor 


The most general renormalizable Lagrangian for scalar mediation consistent with the 
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above assumptions is given by 


CsH = - mlS^S - ^{S^Sf 

+\d^(t)d^^(t) - - ^1^03 _ ^04 

+iqpq — niqqq 

-giuisS'^Scl) - - /iigg0 - i/i2g7^g0, (19) 


where we have suppressed all the SM quark interactions. Similarly, the Lagrangian for vector 
mediation (up to gauge hxing terms) is 

CSG, = S,S'S'‘S - m|S'S - 

+iqlpq — niqqq 

- tg,{S^d^S - d^S^S)G^^ 

-h3{ql^^q)G^ - /i4(g7M7®g)G^- (20) 


2. Spin-\ Dark Matter 

If the WIMP has spin-1 (denoted by y below), then, as in the scalar WIMP case, me¬ 
diation will only occur via scalar or vector mediators. The most general renormalizable 
interactions for the scalar (0) and vector mediator {G^) cases respectively are given below. 


GxH = ixJPx - m^XX 


i3 _ h'2 ,4 

3 ^ 4 ^ 


+iqpq — niqqq 

-Ai0xx - fA20 x7^x - hi4>qq - *^ 20 ^ 7 ^ 9 , 


( 21 ) 
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l^xGq = - m^XX 

+iqlpq — rriqqq 
-AsXT^xG;. - \iXl^l^xG^l 

-h3qX^lqG^ - h^q-ff.'y^qG^. ( 22 ) 

3. Spin-1 Dark Matter 

If the WIMP is a massive spin-1 particle, uncharged mediation to the quark sector can 
occur via a heavy scalar or a vector particle. For the case of vector mediation, there are 
many possible interactions because the Lorentz indices on the vectors afford a more diverse 
set of terms. The general interaction Lagrangian for the scalar mediation case is 

I G.^2 ^ 2 i2 ^0/^1 i3 ^2 ,4 

+ 2('9/^W - 2^4^^ -3 

-\-iq]pq — rUqqq 

-bimx(l)XlX^ - - hicpqq - (23) 

For the case of vector mediation, there are many possible interactions because the Lorentz 
indices on the vectors afford a more diverse set of terms. The Lagrangian is given by 

^xa, = + mi A* A'" - ^(AjA^" 

-15^.5'“'+lrnlGl-X^G,G‘‘f 

-\-iqIpq — rUgqq 

-^-iGliXlxn - |(G^G^)(XtX,) - [thXld.X^'G'^ 

+bQXld^X,G’^ + bje^,,,{X^^d'^XP)G^ + h.c: 

-hG^qYq - h.G^qYxh (24) 

where, for the Lagrangian to be Hermitian, b^ and 67 are complex (this implies a new source 
of CP violation, which will not be considered further here). 
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A. Charged-mediator Lagrangians 


Here we consider the simplest case of mediators that are charged under both the DM 
internal symmetry group and SM gauge groups. This is motivated by the absence of spin- 
I mediators (s-channel processes) in the previous section. Such a mediator, if neutral, is 
forbidden by simultaneous requirements of gauge invariance and renormalizability. Dark 
Matter models with mediators endowed with charges from both DM and SM side have been 
considered in the literature before [57, 58]. The case of a spin-| mediator carrying SU{3)c 
is also motivated by studies of heavy quark models. This allows unique interactions as we 
show below. In particular they necessitate a direct interaction between quarks and WIMPs 
at the level of the Lagrangian. 

1. Scalar Dark Matter 

Scalar WIMPs with a charged scalar or vector mediator do not lead to any Lorentz 
invariant interactions. This is easy to see since both the scalars (or scalar and vector) and 
the quark are required in the (gauge invariant) interaction, but there is no way to contract 
the spinor indices consistently if the mediating particle is a scalar or vector. Therefore, the 
only possibility is that of a spin-1/2 mediator, Q, which acts like a heavy quark. The general 
renormalizable action is given by 


CsQc = d^S^d^S - mlS^S - XsiS^Sf 
+iQIpQ - mqQQ 
+iqlpq — rriqqq 

-{yiSQq+ y 2 SQ'y^q + h.c.), (25) 


where yi and y 2 are again complex. 

2. Spin-\ Dark Matter 

For a spin-1/2 WIMP, both a charged scalar and charged vector mediator exchange can 
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lead to novel interactions. The charged scalar is denoted by <h and the charged vector by 


Cx^q = ixJPx - rrixXX 
+iqlpq — ruqqq 

-{li^^XQ + k^^Xl^q + h.c.), (26) 

^xVq = ^X'^X - mxXX 

+iql/)q — rUqqq 

-{diXl'^qV^ + d 2 XX^X^qVi + h.c.), (27) 


where li,l 2 ,di and d 2 are complex. 

3. Vector DM 

Here again we only have the case of a spin-^ mediated interaction between vector DM 
and quarks (again scalar and vector charged mediators aren’t possible because they don’t 
lead to Lorentz invariant and renormalizable interactions). The general Lagrangian is given 
by 

CxQ, = + mlxlx>‘ - ^(XlX^Y 

+iQlpQ - mgQQ 
+iqlpq — rUqqq 

-{ys^f^Qj^q + y^X^Qx^X^q + h.c.), (28) 


where y^ and y^ are complex. 

IV. NON-RELATIVISTIC REDUCTION OF SIMPLIFIED MODELS 

After integrating out the heavy mediator we replace quark operators with nucleon oper¬ 
ators (see appendix C), take the non-relativistic limit (see appendix B), and match onto the 
operators given in table I. The results of this calculation are presented in terms of the Cj 
coefficients from [54], described in section II, facilitating a straightforward computation of 
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amplitudes and rates. The q’s are given for each of the WIMP spins in tables II, III and IV. 
With this general framework in place we can now easily hnd the leading order NR operators 
for each distinct WIMP-nucleus interaction. One can imagine a series of minimal scenarios 
in which a combination of two Lagrangian couplings that give rise to a direct detection 
signal is non-zero with all others set to zero, and then proceeding in this manner for the 
entire set. Each of these scenarios is listed with its leading operators in table V and with 
all operators generated in table VI. Note that in the case of a complex coupling constant 
we consider purely real and purely imaginary values as separate cases since they produce a 
distinct set of operators. 

TABLE II. Non-zero q coefficients for a spin—0 WIMP 



Uncharged Mediator 

Charged Mediator 

Cl 

CIO 

hfgi 

2 

777-1 

ft> 

-ih^gi 1 2ig4h^ 
nil m-g mg 

ylyi-yly2 
mqms •'T 

■y2yi-yly2XN 

rngms 


TABLE III. Cj coefficients for a spin-i WIMP 


Uncharged Mediator 


Charged Mediator 


Cl 

C4 

C6 

C7 

C8 

eg 

cio 

Cll 

Cl2 


<Ai 
-^- 
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A2mjv 

2hf A4 


2h^ \4 


hJVl 


fe"Ai 
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\2mj^ 


— (12^2—d\di 

4m^ 4mL 

$ V 


/: 


N I I lll2+l\h I dld2+dldi 

+ I 4iii^ + 


8rn^ 




Ijh+Hjh 
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ci|d2— 

2my 
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//J/p ^2^2 I (^2^2 ^1^1 fTliV N 

V Zlm? dm?. / 


/ ^2^2 ^2^1 I (i^ci2?~d2^i \ A TV 
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2ml 
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2m| 


4mt 
d|d2-|-(i2<^l 
4m|, 




//j /2 ^ 2 ^^ ^2 \ a a N ^^1^2 ^ 2 ^^ ^2 H“C^2^1 \ hflpy a 

^ 2m|^ 4m^ '' ^ 2m|^ ArUy ^ 

•( — I ^2^1~^i^2 ^ A /|^/2~^2^1 

V 4m.?. 4m.?_ / m? 


N 


]( ^ 2 ^^ ^2*^1 dj-^d2 \ 771.jY ^TV 

“T /t-v-.2 T 777,^ ./ T 
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47771” 

$ 


4rri^ 


-f ^ 


. / j ^/2 ^2^1 ^iV 


:A^ 


4d 


As described earlier, we hnd that it is important to consider operators beyond those 
incorporated into the standard spin-independent and spin-dependent formalism, i.e. simple 
models exist in which one would infer an incorrect rate in current experiments by not in¬ 
cluding these effects. Also importantly, not all of the NR operators are actually generated at 
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TABLE IV. Ci coefficients for a spin-1 WIMP 


Uncharged Mediator 


Charged Mediator 


Cl 

C4 

C5 

C6 

C8 

eg 

cio 

Cll 

Cl2 

Ci4 

Cl7 

Cl8 


4Im(67)/i; 






m?. 

<p 

Ke(h7)hi 


Re(b6)/i 


mxrtiN 


Ile{be)h^ mjY 


mx 


Re(fe6)fe^ mjv 




mx 


— I 


■ Ke{b 7 )h^ m 






2Im(b7)h^ 

2Re(b6)fef mx , 2Im(fc7)fe^ 

rrig mx ' rriQ 


bih 


N 


hiV , 


mx 


H.e{b'Y)h^ Yfij^ 


mx 

0 


2Ke{bx)h^ 

mx 


mx 

41m(bQ)h^ 

mx 


mx 

4lm(be)h^ mx 




mx 


ylvi-ylvi rN 
mgmx •’T 

n yjys-ylyi xN 

mgmx 

0 

0 

0 

0 

■ yhi-ylyi XN 

mgmx 

^ ylys-yjyi xN 

mgmx 

2 My*-yly3 xn 
mgmx 

0 

0 

2 My3-yly4 xN 
mgmx 


leading order; for example, the operators O 2 , O 3 , O 13 and O 15 are missing at leading order. 
Note that we only consider renormalizable Lagrangians, higher order non-renormalizable 
operators, which are presumably further suppressed. We have also not considered the case 
of kinetic mixing, which could be used to generate anapole interactions [56], because the 
effective interaction doesn’t arise from one mediator exchange. 

While spin independent interactions are a generic feature of direct couplings to quarks 
in our charged mediator cases, it is sometimes possbile to suppress them. In the scalar (and 
vector) WIMP with charged mediator cases, it is possible to suppress the spin independent 
interaction by ensuring that \yi\ = ||/ 2 |(| 2 / 3 | = \yi\) while keeping their relative phases non¬ 
zero (or tt). While these non-minimal scenarios require some hne tuning we include it for 
completeness and label them 1 / 1 , 1/2 and 2 / 3 , 1 / 4 . 

Aside from scalar WIMPs, each particular spin produces some non-relativistic operators 
that are unique to that spin. Also, importantly, the operators Oi and Oio are generic to 
all spins. In hve cases relativistic operators generate unique non-relativistic operators at 
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leading order. Therefore distinguishing WIMP scenarios in these cases reduces to experi¬ 
mentally discerning between these operators (see also [59]). Given the likely low statistics 
of any detection in upcoming direct detection experiments, sub-leading operators are not 
likely to contribute enough to provide any further discriminating power. 


V. OBSERVABLES 


The principle observable in direct detection experiments is the differential event rate. 
Since the incoming WIMPs originate in the galactic halo, one must average over the WIMP 
velocity distribution, f{v), which we assume for the purposes of this paper to be Maxwell- 
Boltzmann, 


dR 

dEn 


Nt 


PxM 

271771 ^ 



/(^) 


Ptot dv 


(29) 


where we use the value = 0.3GeV/cm^ for the local dark matter density, Nt is the number 
of nuclei in the target and Ptot can be calculated from the amplitude M. in Eq. 14 


Ptot — 


2jy + 1 2jAr -|- 1 


E l-wp 


spins 


(30) 


Thoughout this work we use the mathematica package supplied in [54] to calculate rates. To 
determine the leading order operator which arises from a given relativistic scenario we hrst 
plot the rate for each of the NR operators in xenon-131. To simply compare the operators 
we set the c, coefficients to be the same and normalized the overall rate to that of Oi, 
see Fig. 1. Since operators are either zero, hrst or second order in momentum transfer q 
or velocity P-*-, the relative strengths of the operators span 16 orders of magnitude. This 
is an important point to keep in mind when hnding the leading operator, as sometimes a 
term which appears to be higher order in q can dominate the non-relativistic reduction. For 
example in the scenario, one hnds that dominates over the Oq and Ou which 

contain powers of q within the operators. 


Since the Lagrangians we have considered are not tied to specihc complete and consistent 
particle physics models, the mediator masses are not hxed in advance and thus specihc 
event rates are not predicted in advance. Glearly one requires a rate that is low enough 
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FIG. 1. The relative strength of event rates for a 50GeV spin-g WIMP in xenon for each of the 
non-relativistic operators in table I, where the coefficients of each operator are set to be equal 

to evade the current experimental constraints. For example, a 50 GeV WIMP producing 
10 events per tonne per year is sufficiently low to evade the bounds from LUX [21]. For 
demonstration purposes we set the couplings to 0.1 (or O.lf for imaginary) in the various 
Lagrangians and hnd a mediator mass that will produce 10 events/t/y in the signal region 
for xenon (5 — 45keV). The calculated masses are given in table V. It is perhaps telling that 
the mediator masses span 6 orders of magnitude, from just a few GeV up to a PeV. While 
it is unlikely that a full model of thermal relic dark matter could be built around all of 
these Lagrangians, it is nevertheless a useful metric to estimate the relative strength of the 
different nuclear responses to each of the operators. 

In Figs. 2, 3, 4 and 5 we have plotted rates for two common targets. For simplicity and 
again for demonstration purposes, we only plot the rates for a single isotope of both ger¬ 
manium and xenon. The choice of isotopes, ^^Ge and ^^^Xe, was made to ensure sensitivity 
to spin-dependent responses. As can be seen in the hgures, many operators produce rates 
with similar recoil energy dependence in the same target, but different nuclei can have very 
different responses to the various operators [53]. Thus a complementary choice of nuclear 
targets can provide important discriminating information. 

To illustrate this discriminating power we plot the ratio of the rates in xenon and ger¬ 
manium in Fig. 5 and 6. We choose to only present ratios for the uncharged mediator 
cases of spinor and vector WIMPs since the other cases produce trival results (all operators 
being spin independent). To estimate the effect astrophysical uncertainties will have on 
discriminating between operators, we plot the rate for a range of astrophysical parame- 
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ters from vq = 200m/s, and Vesc = 500m/s (lower) to no = 240m/s and Vesc = 600m/s 
(npper). The nncertainty in the dark matter density does not appear since we are con¬ 
sidering the ratio of rates. Given the vastly different energy dependence of the ratio of 
rates of each scenario the astrophysical errors do not completely inhibit their identification. 
Furthermore, operators Og and Ou, produced in scenarios and h 4 bf^ respectively, 

remain indistingushable when considering the ratio of rates. While it appears that in prin¬ 
ciple almost every operator is discernible, in practice isotopically impure targets and low 
statistics will further complicate the situation and provide limits on practical discrimination. 




- yi 

- 72 

- 7172 


FIG. 2. Rates for a 50GeV spin-0 WIMP in xenon (solid) and germanium (dashed) with uncharged 
(left) and charged mediators (right), assuming mediator mass of ITeV and 0(1) coupling constants. 




FIG. 3. Rates for a 50GeV spin-^ WIMP in xenon (solid) and germanium (dashed) with uncharged 
(left) and charged mediators (right), assuming mediator mass of ITeV and 0(1} coupling constants. 
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FIG. 4. Rates for a 50GeV spin-1 WIMP in xenon (solid) and germanium (dashed) with uncharged 
(left) and charged mediators (right), assuming mediator mass of ITeV and 0(1) coupling constants. 




FIG. 5. Rates (left) for a 50GeV spin-1 WIMP in xenon (solid) and germanium (dashed) with 
uncharged mediators and imaginary couplings, assuming mediator mass of ITeV and 0{1) coupling 
constants. Also shown is the ratio of rates in xenon and germanium (right). 

VI. CONCLUSION 

The analysis we have given here bnilds on previous analyses to provide, in generality, a 
roadmap to use event rates in direct dark matter detectors to constrain fundamental dark 
matter models. We have outlined the steps needed to go from fundamental Lagrangians, hrst 
to relativistic operators, then to non-relativistic operators, and finally to produce nuclear 
matrix elements. In the process several significant facts have been elaborated. 

• Not all possible non-relativistic operators contributing to nuclear matrix elements in 
direct detection will arise from simple UV complete dark matter models. 

• Aside from scalar WIMPs each particular spin produces some non-relativistic operators 
that are unique to that spin. 
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— /7ibi 
- /72bi 

— llAbs 

— 

— hAbe^^ 

— h^t>7^^ 


FIG. 6. Ratio of rates in xenon and germanium, illustrating the discriminating power of having 
multiple nuclear targets. For a 50GeV spin-| WIMP with uncharged mediator (left) and a 50GeV 
spin-1 WIMP with uncharged mediator (right), the shaded regions show the upper and lower 
bounds due to the astrophysical parameters 

• Two non-relativistic operators, Oi and (9io, are ubiquitous and arise for all WIMP 
spins we have explored. 

• In 5 scenarios, relativistic operators generate unique non-relativistic operators at lead¬ 
ing order. 

• Two new non-relativistic operators not previously considered within the context of 
the full array of allowed nuclear responses arise at low energies if spin-1 WIMP dark 
matter is allowed for. 

• While the different operators that can contribute to event rates in detectors using 
specihc elements or isotopes cannot be distinguished on the basis of their impact on the 
differential event rates in these detectors, they can produce radically different energy 
dependence for scattering off different nuclear targets. Thus, a complementary use 
of different target materials will be necessary to reliably distinguish between different 
particle physics model possibilities for WIMP dark matter. 

While current detectors have only yielded upper limits, with new generations of larger 
detectors with greater energy resolution and lower thresholds coming online, the search for 
WIMP dark matter has never been so vibrant and promising. The tools we have provided 
here should help experimenters to probe the most useful parameter space, to interpret any 
non-zero signals in terms of constraints on fundamental models, and should allow theorists 
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who build fundamental models to frame predictions in an accurate and simple way so that 
they might be directly compared with experiment. 
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TABLE V. Leading order operators which can arise from the relativistic Lagrangians considered in 
this work, the column ‘T terms’ gives the non-zero couplings for that scenario. Each row represents 
a possible leading order direct detection signal. A ‘f’ indicates that the mediator is charged. The 
’Eqv. Mm column gives the mediator mass required for each scenario to produce ~10 events 
t~^yr~^keV~^ in xenon, with couplings set to 0 . 1 . 


WIMP spin 

Mediator spin 

C terms 

leading NR operator 

Eqv. Mm 

0 

0 


hi,gi 

Oi 

13 TeV 

0 

0 


h2,gi 

Oio 

14 GeV 

0 

1 


h4,gA 

Oio 

8 GeV 

0 

it 

2 


yi 

Oi 

3.2 PeV 

0 

it 

2 


2/2 

Oi 

3.2 PeV 

0 

it 

2 


yi,y2 

Oio 

41 GeV 

1 

2 

0 


hi,Xi 

Oi 

12.7 TeV 

1 

2 

0 


h2, Ai 

Oio 

293 GeV 

1 

2 

0 


hi,X2 

On 

14 GeV 

1 

2 

0 


^2, X 2 

Oe 

1.9 GeV 

1 

2 

1 


hs, A 3 

Oi 

6.3 TeV 

1 

2 

1 


/ 14 , A 3 

O 9 

6.4 GeV 

1 

2 

1 


/ 13 , A 4 

Os 

180 GeV 

1 

2 

1 


/ 14 , A 4 

O 4 

135 GeV 

1 

2 

Ot 


h 

Oi 

7.1 TeV 

1 

2 

ot 


I 2 

Oi 

5.5 TeV 

1 

2 

it 


di 

Oi 

5.9 TeV 

1 

2 

it 


d2 

Oi 

6.7 TeV 

1 

0 


hi,bi 

Oi 

13 TeV 

1 

0 


h2,bi 

Oio 

10 GeV 

1 

1 


hi, 65 

Oio 

5.1 GeV 

1 

1 

/13 


05(017) 

5.5 GeV(23 GeV) 

1 

1 

/14 

,bf%b]r) 

09 ( 0 i 8 ) 

3 GeV(4 .6 GeV) 

1 

1 


,bf%b^^) 

Oii( 08 ) 

186 GeV(228 GeV) 

1 

1 

/14 

,bf%b\^) 

04(04) 

78 MeV (172 GeV) 

1 

it 

2 


2/3 

Oi 

3.2 PeV 

1 

it 

2 


2/4 

Oi 

3.2 PeV 

1 

it 

2 


2/3, 2/4 

On 

120 TeV 
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TABLE VI. List of scenarios with leading operators colored by which are distinguishable via the 
ratio 




O 2 

O 3 

O 4 

q^Oi 

Oe 

Oe 

O 7 

Os 

O 9 

Oio 

Oil 

O 12 

Ol3 

Ol4 
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Ol7 
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/ 
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/ 

(y4) 
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/ 





/ 
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/ 

/ 
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* indicates the purely imaginary scenario for that coupling 
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Appendix A: Vector Dark Matter 


If the WIMP has spin 1, we hnd two extra operators that haven’t been considered pre¬ 
viously. Specihcally, the operators depend on the symmetric combination of polarization 
vectors, Sij = | (ejej + This necessitates a modihcation to the WIMP response func¬ 

tions by hrst modifying the £ coefficients given in Eq. 13. Based on our non-relativistic 
reduction for vector dark matter, the Lagrangian for vector dark matter and the nucleus, 
interacting via an uncharged scalar or vector mediator can be written in general as: 


£^vector — CiOi + 04^04 + C5O5 + CgOg + CgOg + CigOiQ + CuOn + C14O14 + C17O17 -j- CigOis 

(Al) 


where we’ve dehned O 17 = ^ ■ S ■ v± and Ois = ■ S ■ and the Cj’s are given in 

table IV. To decompose these new operators we replace v-^ with the target velocity and the 
internucleon velocities and sum over nucleons. O 17 can then be put into the form 

A 


iriN 

Ois can be expanded as 


Vrpe 


— iq.x 


2 = 1 


i _ 'W -^7 -I- 

Z * ^ t y^nt 


(A2) 


ZrriN 


(A3) 


Together, all the terms of elector give rise to the following £ factors from Eq. 13, 


= c[ -1-z 


mN 


X 4 ■ S^cl + {v^ ■ S^)cl + i 


mN 


■ r^n + * 


mN 


■S-vl\c\j 


lo^ = 


11 = 0 


q 




(A4) 


Im — i 


mN 


xS^]cl-S^cl-z 


■5 c 


mN 


17 


^5 - i 


1 /. <? 


X Co + - U 


mN ® ' 2 


rriN j 2^"^ i 2mAr 


• c[4 + 


1 /. <? 


■5 c 


mN 
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Based on the £’s above, the coefficients of the various nuclear responses are found by squaring 
the amplitude and then summing over spins. To simplify calculations, we choose a convenient 
basis for polarization vectors, e| = 5®. Recall that the spin can then be written as the anti¬ 
symmetric combination iSk = eijke\ej. The WIMP responses unique to the vector case are 
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then given by: 


prr _ T T 


Rl^„ =0 

M ^ 0 



^ ±2 r r' , ±2 r r' , 

-^Vt C5C5 + Vt CgCg + — 


„2„,±2 ^ 

, , , Q T t' 

-CiiCii + ^^2 Ci7C;^7 


Af 


4m 


Af 


=0 


-RhT/ — — C4C4 + 
0 


PTT _ _„T T 

n^, gC4C4 


<f r t' , (f T r' 




6m 


■CqCq + 


2 I 2 

q 


^ ' 2m2 


C[4C[4 + 


DTt' _ 

— ~ 




N 

o 

9 


24m^ 


'^18*^18 


m 


= i ( 


N 
5'^4 


+ g^2 C17C17 


N 


T T T T 

- CgCg 


(AS) 


Appendix B: Non-relativistic Reduction 


We find effective relativistic interaction Lagrangians by integrating out heavy mediators. 
We only keep the leading order interactions (suppressed by m or m?). To the right of each 
operator is their non-relativistic reduction expressed in terms of the operators in table I with 
the coefficient derived from the Lagrangian parameters along with the relevant nucleon form 
factor. As multiple operators can have the same non-relativistic limit, it is important to 
include the nucleon form factor at the relativistic level. If this is not performed, erroneous 
cancellations can occur. 

For free spinors we use the Bjorken and Drell normalization and 7 matrix conventions. 
In the non-relativistic limit we make the following replacements: 

I 5 






X ^ 




^/mx 


E -I- m. 


2m^ 



(Bl) 


where s = 1, 2, 3 are the different polarization states of the vector. .^ = (1 0)^ is the left 
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handed Weyl spinor. The following Fierz transformation and gamma matrix identites were 
usefnl in the charged mediator cases, (a sign difference was fonnd in the hnal identity when 
compared with [60]): 


{qx){x<i) =- 
{qi^x){xi^(i) =- 
{qx){x'i^(i) =- 

{Q'l^lX){xY(l) =- 

=- 

{Ql^lX){X'l^'l^(l) =- 


1 

'4 

1 

'4 

1 r 

■4 [' 


mxx + qYqxi^x + ^Q(^'"'QX(^^,uX - Ql'"YQXl^,Yx + qYqxYx 




1 


qqxx + qi^qxi^x - qi^qxi^x + qi^i^qxivn^x + ^q(^^"qx(^^yX 
qqxi^x + qi^qxx - ql^qxl^ll^x + qi^i^qxii^x + i(^t,uaf^q(^^’'qw°‘^x 


qqxx - qi^qxi^x - ^qi^qxii^x - ^qi^x^qxifix^x 


-qqxx + qx qxi x - ^qx^qxxt^x - ^qx^x qxxt^x x 


qqxx x-qx qxx + ^qx^qxx^n x + -jp^x qxxi^x 


(B2) 


a- 


- 7 = = 2^' 


pupa 


a. 


pa 


(B3) 


All of the following operators are collected in terms of the coefficients of the NR operators, 
Q, in tables II,III and IV. 


TABLE VII. Non-relativistic reduction of operators for a spin-0 WIMP 
Scalar Mediator 

{S^S){qq) Oi 

(gt 5 )(g^ 5 g)(^) 010 

Vector Mediator 


i{S^d,S-d^S^S){q^^^^S) 010 


Charged Spinor Mediator 
{S^S){qq) 

{S^Sm^q) 


ylyi-yly2 ^nq 

rriQms ^ 


^ yW-yW -^NQ 

rriQms 

























TABLE VIII. Operators for a spin-^ WIMP via a neutral mediator 
Scalar Mediator 



Vector Mediator 

Xl^Xqi.l^q (-^) (-^7 + ^09) 

Xl^fxqi^q (-^) (^8 + O9) 

xi^i^xqiiii^q —^ 
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TABLE IX. Non-relativistic reduction of operators for a spin -2 WIMP via a charged mediator 
(after using Fierz identities) 

Charged Scalar Mediator 


XXQQ 


xxqx^q 


XTxqq 


xi^xqi^q 


xTxql^lq 

xi^i^xqiiiq 

xi^xqxt^'i^q 

xi^i^xqiiii^q 


xcr^''xq(^liiyq 


ZjL fN /'o 
Ami ^ 


^ ZjZl z(z2 mjv fN /O . 
^—^zps—TA^JTq^n 


4mT 

•i* 


Aq^Oe 


l-^l\ ^2^2 '^N 
4m| 

%|^A/7(08 + O 9 ) 

zlZ2+4h 

2m? 


- ^ 09 ) 


_ 4h+4L ^^0 

9 

4L-4h 

mt. 9 
$ 

- i^Oii + 4 O 12 ) 


Charged Vector Mediator 


xxqq 

xxqi^q 


XI xqq 
xn^xqi^q 


xi^xqitiq 

xn^i^xqiiiq 

xi'^xqxt^x^q 

xn^i^xqiiii^q 


C?2^2 (l^d\ pJXJ 

4m ^ 

A '-AT/z-j 


4m 


V 


.c/2(^l cl^d2 ITljsi -pN 


4mi 


V 

^ 2 ^ 2 —m^v /\ zN 
4m^ 

d^fi2+£^Jc?i \rNrn 

+ o,) 

a!^d2+cZ(iil \Nrn 
2m^ ^9 ‘-^4 
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TABLE X. Non-relativistic reduction of operators for a spin-1 WIMP 


Scalar Mediator 
XlXf^qq 
XlXf^qj^q 
Vector Mediator 



b,h^\ 

) 

bih^ 


Oi 


Ow 


{Xld^X’' - d^XlX^iq^f^q) 

{Xtd^X'^ - 

a,(A"tx^ + AtA")(^ 7 '^g) 

a,(A"tx^ + AtA")(^7V'i) 

d.{X^^X^-Xlxn{qi^q) 

a,(A"tx^-AtA")(g7V9) 

e^^upa [x'^^dPX- + {q^^^q) 

epupa [x'^^dPX'^ + X’'dPX^^) {q^^^j^q) 
epypa - X'^dPX^^) {q^Pq) 

epupa (x^'^dPX- - X'^dPX-^) {q-fP-f^q) 


—Sbsh^ mjy 

mx 


f Re{b6)h^ rax 
y mx 

2 Ke{be)h^ mx 


Ow 

(O 5 + 06 - 
09 


o. 


mx 


4Im(b6)fef mjv 


mx 


( 41ra{bQ)h^ mx 

[ 

^Re(fe 7 )/i^ 


Ol7 




mx 


'3 mx 


mx 


O 18 

On 


/ Re{b7)h^ mx 
mx 

( 21 m(b'Y)h^ 


{ 


) (’Sij+04-iSf06 - 20l4) 


/ 4Im(fe7)fei' 


(Os + O 9 ) 
O 4 


Charged Spinor Mediator 

iXlX,){qrYq) 

iXlX,){qjPY^^q) 


yjya- 


■vlvA 


mgmx 


(/f^Oi+25^04 


y\y3- 


■ylyi 


(iAf Oio + iV^Oii - 2iV"Oi2 - 2iV"Oi8) 


7N 


mQmx 


Nt 


7 N 7 


Appendix C: Quarks to Nucleons 

To go from the fundamental interactions of WIMPs with quarks to scattering from point- 
like nucleons, one must evaluate the quark (parton) bilinears in the nucleons. For a full 
discussion see the appendix of [60] and [61]. We write the nucleon couplings in terms of the 
quark couplings times a form factor (in the limit of zero momentum transfer): The scalar 
bilinear for light quarks can be evaluated from 

{N\ mgqq\N) = mjv/Tg (Cl) 
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{No 

mqqq\Ni) - 

f^qNN 

{No 

qi^q\Ni) - 


{No 

qj^qlNi) - 


{No 

qi^l^qm - 


{No 

qa^^'^q\Ni) - 

Syaf^’^N 


while for the heavy quarks 


{N\ m^qq\N) = ( 1 - frq 

' q=u,d,s 


(C2) 


Summing over all the quarks one finds 


h^= y — y 

q=u,d,s y q=c,b,t h 


(C3) 


The psuedo-scalar bilinear was recently revisited in [61]: 


Ar= E l4Af-AG^ E — 

^ '^a 

q=u,d,s q=c,bd h 


(C4) 


The vector bilinear essentially gives the number operator: 


y = 


2/13 + ^3 N = p 
hn + 2hi N = u 


(C5) 


The psuedo-vector bilinear counts the contributions of spin to the nucleon (note that 
sometimes this coupling has a, Gp factored out to make it dimensionless) 


y = y hiA, 


(C 6 ) 


q=u^d^s 


Throughout this paper the following values are used (it should be noted that there are 
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large uncertainties in these values) [60, 61]: 


£71 

JTu 

= 0.014 

II 

-- 0.02 

£71 

JTd 

= 0.036 

II 

= 0.026 

£71 

JTs 

= 0.118 

fp _ 
JTs - 

: 0.118 

<1 

= -0.427 

K = 

0.842 

<1 

= 0.842 

A^ = 

-0.427 

a : 

= -0.085 

A? = 

-0.085 

Ah" 

= - 108.03 


= 110.55 

AJ" 

= 108.60 

A(F -- 

= -107.17 

<1 

= -0.57 

AsP -- 

= -3.37 

AG" 

= 35.7MeV 

AGP 

= 395.2MeV 


Assuming a universal coupling of the mediators to all quarks, the nucleon level couplings 
can then be written as. 


h^=f¥h 
=A^h2 
hf =M^h^ 
hf =A^h4 

where we have dehned, 

= 11.93 = 12.31 


A" = 

-0.07 

A^ = -0.28 

AA" = 

3 

AfP = 3 

A" = 

0.33 

AP = 0.33 

5" = 

0.564 

6 P = 0.564 


(C8) 


(C9) 
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This introduces a small amount of isospin violation, and it is known that relaxing the assump¬ 
tion of universal couplings to quarks can lead to interesting isospin violating effects [61, 62], 
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